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Abstract 

To  support  the  impedance  analysis  in  a  disk- type  SOFC  under  power  generation,  a  numerical  simulation  procedure,  whose  simulation 
outputs  the  total  impedance  of  an  actual-size  disk- type  SOFC  cell  for  the  input  of  the  local  impedance  of  cell,  has  been  proposed.  In 
the  simulation,  the  quasi-one-dimensional  radial  flow  and  mass  transport  problem  with  cell  reactions  is  solved.  Transient  distributions 
of  electromotive  force,  current  and  concentrations  of  chemical  species  in  the  cell  are  calculated  to  obtain  the  total  cell  impedance  with 
parameters,  such  as  local  impedance  of  the  cell,  ac  frequency,  fuel  utilization,  gas  flow  rate,  etc.  A  parametric  simulation  study  for  the 
simplest  case  under  typical  experimental  conditions  has  been  carried  out.  It  is  found  that  the  Cole-Cole  plot  for  the  total  cell  impedance 
shows  a  capacitive  semicircle,  even  if  a  pure  resistance  is  assumed  as  the  local  impedance  of  the  cell.  The  diameter  of  the  impedance 
semicircle  changes  with  fuel  utilization  and  supplied  gas  flow  rate,  being  affected  by  the  diffusion  in  the  flow  direction.  The  characteristic 
frequency  of  the  impedance  semicircle  depends  on  the  transit  time  of  the  gas  through  the  cell  which  is  also  affected  by  the  diffusion  in  the 
flow  direction.  It  is  also  suggested  that  the  total  cell  impedance  appears  as  a  convolution  of  this  capacitive  semicircle  impedance  and  local 
impedance  of  the  cell. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  ac  impedance  analysis  is  one  of  the  most  effective 
methods  not  only  for  characterization  of  materials  and  elec¬ 
trodes  in  electrochemical  systems,  but  also  for  obtaining 
characteristic  data  necessary  for  system  design  and  simula¬ 
tions.  The  equivalent  circuit  obtained  from  the  impedance 
analysis  would  be  especially  useful  for  the  simulation  of 
the  dynamic  performance  of  such  systems.  Various  types 
of  solid  oxide  fuel  cells  have  been  under  development  in 
the  USA,  Europe  and  Japan,  and  the  concept  is  already  in 
the  demonstration  stage  for  commercialization  at  Siemens 
Westinghouse  Power  Corporation  [1],  Sulzer  Hexis,  Ltd. 
[2],  etc.  However,  impedance  analyses  of  commercial  size 
SOFCs  under  practical  power  generation  conditions  have 
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been  scarcely  reported  [3,4].  The  main  reason  for  this  is 
that  the  measurement  and  the  interpretation  of  impedance 
in  a  full-scale  cell  is  rather  difficult,  due  to  the  nonuniform 
current  distribution  in  the  cell  and  the  unknown  correlation 
between  the  local  impedance  of  the  cell  and  the  total  cell 
impedance. 

Disk-type  SOFCs  have  currently  been  developed  by 
Sulzer  Hexis,  Ltd.  [2]  and  Mitsubishi  Materials  Corporation 
in  collaboration  with  the  Kansai  Electric  Power  Company, 
Inc.  [5,6].  In  a  disk-type  axially  symmetric  cell  (e.g.  Mit¬ 
subishi  Materials  type),  the  impedance  measurement  and 
analysis  may  be  greatly  simplified  [6] . 

Many  experimental  investigations  on  impedance  of  the 
SOFC  cell  have  been  conducted  using  small-size  cells 
[2,7-12].  Generally,  the  impedance  of  the  SOFC  cell  is 
related  with  many  factors,  such  as  ionic  and  electronic  con¬ 
ductivities,  diffusivities  of  the  chemical  species  in  the  elec¬ 
trolyte  and  electrodes,  reaction  activation  over-potentials 
at  the  anode  and  cathode,  and  diffusion  of  the  chemical 
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Nomenclature 


c 


Cd 

cg 

cwg 

D 


TAell 

^emf 


F 

AG 


o 

h2o 


h 

II 

Jr 

N 

P 

r 

R 

T^cell 

Ri 

Rr 

Rs 

Ry/g 

rwg 

?d 

t 

T 

Uf 

ue  ff 

v 

V 

Vac 

Vdc 

Vcell 


molar  density  (mol/cm3) 
double  layer  capacity  per  unit  area  (F/cm2) 
parallel  capacity  for  gas  conversion 
impedance  (F) 

normalized  parallel  capacity  for  whole-cell 
gas-conversion  impedance  (F/cm2) 
diffusion  coefficient  (cm2/s) 
cell  diameter  (cm) 
electromotive  force  (V) 

Faraday  constant  (C/mol) 

Gibbs  free  energy  change  in  the  reaction 

H2  +  ^C>2  ^  H2O  at  the  standard  pressure 
(0.1  MPa)  and  reference  state  (J/mol) 
channel  effective  height  (cm) 
load  current  (A) 

current  density  for  reaction  (A/cm2) 

volumetric  flow  rate  (cm3/s) 

total  pressure  or  partial  pressure  (0.1  MPa) 

radius  (cm) 

gas  constant  (J/mol/K) 

local  cell  resistance  per  unit  area  (Q  cm2) 

electronic  resistance  per  unit  area  (Q  cm2) 

ionic  resistance  per  unit  area  (Q  cm2) 

parallel  resistance  for  gas  conversion 

impedance  (Q) 

reaction  resistance  per  unit  area  (Q  cm2) 
serial  Ohmic  resistance  per  unit  area  (Q  cm2) 
normalized  parallel  resistance  for  whole-cell 
gas-conversion  impedance  (Q  cm2) 
time  constant  for  whole-cell  gas-conversion 
impedance  (s) 

time  constant  for  diffusion  in  z-direction, 
D/h2  (s) 
time  (s) 

temperature  (K) 

fuel  utilization 

effective  fuel  utilization 

molar  average  flow  velocity  (cm/s) 

controlled  volume  by  electrode  reaction 

(cm3) 

amplitude  of  ac  voltage  (V) 

dc  voltage  (V) 

cell  terminal  voltage  (V) 


Subscripts 
a  air 

f  fuel 

H2  hydrogen 

H2O  water  vapor 

M  Ff>,  H2O  or  O2 

O2  oxygen 

r  radial  axis  direction 

z  z-axis  direction 


species  in  the  gas  flow  channel.  It  also  depends  on  the 
cell  construction  and  working  conditions.  If  the  basic  cell 
design,  such  as  the  sizes,  materials  and  structures  of  the 
electrolyte  and  electrodes  and  the  morphologies  of  their 
interfaces  is  given,  the  impedance  may  be  a  function  of  the 
working  conditions,  such  as  temperature,  pressure  and  the 
compositions  of  the  fuel  and  oxidizing  agent. 

This  study  deals  with  a  simulation  method  which  links 
the  local  impedance  of  the  cell  with  the  total  cell  impedance 
at  the  power  generation  under  the  assumption  that  the  lo¬ 
cal  impedance  is  given  by  a  function  of  the  frequency  with 
some  parameters,  such  as  temperature,  total  pressure  and  the 
compositions  of  fuel  and  oxidizing  agent.  In  this  study,  a 
numerical  simulation  procedure  has  been  proposed  in  order 
to  support  the  impedance  analysis  in  a  disk-type  commer¬ 
cial  size  SOFC  under  power  generation.  In  the  simulation,  a 
quasi-one-dimensional  radial  flow  and  mass  transport  prob¬ 
lem  including  cell  reactions  was  numerically  solved  as  a 
function  of  voltage  and  direct  current  superimposed  with  a 
small  alternating  current.  With  a  given  function  of  the  lo¬ 
cal  impedance,  transient  distributions  of  the  current,  electro¬ 
motive  force  and  concentrations  of  chemical  species  in  the 
cell  were  calculated,  and  the  total  output  current  as  a  func¬ 
tion  of  time  was  obtained.  Thus,  the  total  cell  impedance 
was  calculated  by  using  the  Fourier  transforms.  A  paramet¬ 
ric  simulation  study  for  the  simplest  case  under  typical  ex¬ 
perimental  conditions  was  also  carried  out  by  changing  pa¬ 
rameters,  such  as  applied  ac  frequency,  load  currents  (fuel 
utilization),  cell  dimensions,  supplied  gas  flow  rates  and  the 
local  impedance  of  the  cell,  and  the  results  are  discussed. 


2.  Simulation  model 

2.7.  Model  assumptions  and  discussions 

The  schematic  diagram  of  a  disk-type  SOFC  as  a  sim¬ 
ulation  model  is  shown  in  Fig.  1.  The  single  cell  consists 
of  a  composite  disk  made  of  a  solid  electrolyte  with  elec¬ 
trodes  on  both  sides,  structured  current  collectors  and  metal 
separators  for  positive  and  negative  electrodes.  Hydrogen  as 
fuel  and  air  as  oxidant  are  supplied  through  orifices  located 
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Fig.  1.  Schematic  diagram  of  a  disk-type  SOFC. 
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at  the  center  of  each  separator,  and  the  structured  current 
collectors  aided  the  uniformly  distributed  axially  symmetric 
gas  flow.  It  is  assumed  that  the  solid  electrolyte  and  elec¬ 
trodes  and  the  current  collectors  are  homogeneous  with  uni¬ 
form  thicknesses.  Although  cell  reactions  take  place  on  the 
electrolyte  surface,  it  is  assumed  that  the  concentrations  of 
the  chemical  species  are  constant  along  the  z  direction  due 
to  perfect  mixing  by  diffusion.  This  assumption  is  valid  for 
the  transient  behaviors  with  the  characteristic  time  constant 
larger  than  the  diffusion  time  constant  Td  =  D/  h 2,  if  the 
parameter  =  vY0h2/r0D  which  is  a  rough  ratio  of  the  ra¬ 
dial  flow  convection  term  and  the  ^-direction  diffusion  term 
in  the  molar  transport  equation  is  sufficiently  small,  where 
vY0  is  the  radial  flow  velocity  at  the  cell  outlet,  r0  the  cell 
outer  radius,  h  the  height  of  the  gas  flow  channel,  and  D  the 
diffusion  coefficient.  The  precise  discussion  about  this  pa¬ 
rameter  is  given  in  the  appendix.  It  is  also  assumed  that  the 
pressure  drop  in  the  channel  is  small.  Therefore,  the  total 
pressure  is  constant.  For  further  simplicity,  the  temperature 
is  assumed  to  be  uniform  in  the  cell  in  this  paper,  though  it 
should  be  determined  by  solving  the  heat  transfer  problem 
in  the  cell  in  a  future  work.  This  assumption  can  be  approx¬ 
imately  realized  by  externally  heating  the  separators  during 
the  experiment. 

The  conductivity  of  the  current  collector  is  sufficiently 
large,  hence,  the  cell  potential,  Vceii  between  the  cathode  and 
anode,  is  constant  everywhere  in  the  collectors.  The  thick¬ 
nesses  of  the  electrodes  and  electrolyte  are  sufficiently  thin 
relative  to  their  areas,  and  the  currents  in  the  r-direction  are 
much  smaller  than  the  total  current.  Therefore,  it  is  assumed 
that  the  reaction  current  density  and  the  output  current  den¬ 
sity  are  locally  determined  by  the  local  impedance  with  the 
local  electromotive  force  and  the  cell  voltage. 


(a)  (b) 


Fig.  2.  Equivalent  circuits  for  local  cell  in  SOFC:  (a)  simplest  circuit;  (b) 
considering  reaction  resistance  with  double  layer  and  electronic  conduction 
in  solid  electrolyte. 
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where  c  is  the  molar  density,  v  the  molar  average  flow  ve¬ 
locity,  D  the  diffusion  coefficient,  r  the  radius,  JY  the  cell 
reaction  current  density,  h  the  effective  channel  height.  Sub¬ 
scripts:  H2,  H2O,  O2,  f,  a  the  hydrogen,  water,  oxygen,  fuel, 
air,  respectively.  The  cell  reaction  current  density  JY  is  lo¬ 
cally  calculated  by  EQ mf  and  Vceii  from  the  given  function 
of  the  local  impedance  as  well  as  the  output  current  density 
Jo,  as  previously  mentioned.  For  the  simplest  example,  the 
local  impedance  consists  of  a  pure  resistance  Rce\\  as  shown 
in  Fig.  2a,  JY  and  Jo  are  given  by 


T^emf 


Fcell 


^cell 


2.2.  Governing  equations 


Basic  equations  governing  the  system  are  as  follows: 
The  electromotive  force  EQ mf  can  be  obtained  from 


T^emf 


^^h2o  RT.  PihPol 

- — | - In - - 

IF  IF  pu2o 


where  AG^0  is  the  Gibbs  free  energy  change  in  the  reac¬ 
tion  H2  +  5O2  ^  H2O  at  the  standard  pressure  and  refer¬ 
ence  state,  F  the  Faraday’s  constant,  R  the  gas  constant,  T 
the  temperature,  pn2  the  hydrogen  partial  pressure,  po2  the 
oxygen  partial  pressure,  pu2o  the  water  partial  pressure. 
The  cell  voltage  Vceii  is  given  by 

Vceii  —  Vdc  +  V^c  sin(27 zft)  (2) 


where  Vdc  is  the  dc  cell  voltage,  Vac  the  amplitude  of  applied 
ac  voltage, /the  frequency,  t  the  time. 

Transport  equations  for  the  molar  densities  of  hydrogen, 
water  and  oxygen  considering  the  diffusion  in  the  flow  di¬ 
rection  are 


Furthermore,  the  continuity  may  be  described  as 


- rv  f  =  0 

r  dr 

Id  RT  JY 

- rna  =  - 

r  dr  p  4FhSL 


(7) 

(8) 


where  p  is  the  total  pressure.  Finally,  the  relationship  be¬ 
tween  the  partial  pressure  pm  and  the  molar  density  cm  of 
the  chemical  species  M  is  given  as  follows: 


Pm  =  cmRT 


2.3.  Local  impedance  of  cell 

The  function  of  the  local  impedance  is  easily  obtained, 
if  the  impedance  can  be  expressed  by  an  equivalent  cir¬ 
cuit.  According  to  the  impedance  analysis  of  SOFC,  the 
impedance  changes  with  frequency  between  several  MHz 
and  sub  mHz,  and  it  has  been  expressed  by  the  equivalent 
circuits  of  the  combination  of  resistive  elements  and  capac¬ 
itive  elements,  inductivity  elements,  etc.  [6-12].  In  this  pa¬ 
per,  the  impedance  of  the  frequency  range  below  100  Hz  is 
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monitored,  where  the  impedance  may  be  affected  by  the  gas 
flow. 

To  examine  the  effect  of  the  gas  flow,  the  simplest  equiva¬ 
lent  circuit  consisting  of  a  constant  pure  resistance  as  shown 
in  Fig.  2a  is  used  as  the  local  impedance.  Furthermore,  the 
solid  electrolyte  has  electronic  conductance  as  well  as  ionic 
conductance.  Thus,  in  order  to  examine  the  effect  of  the 
mixed  conductivity  in  the  electrolyte  as  well  as  the  com¬ 
plex  local  impedance,  an  equivalent  circuit  as  shown  in 
Fig.  2b,  which  includes  a  representative  RC  parallel  cir¬ 
cuit  with  reaction  resistance  Rr  and  double  layer  capaci¬ 
tance  Cd,  electronic  resistance  Re,  ionic  resistance  R{  and 
series  ohmic  resistance  Rs  is  also  used  for  comparison  and 
discussion. 


2.4.  Numerical  method  and  conditions  in  calculation 


The  previously  described  partial  differential  equations 
were  numerically  solved  by  employing  the  implicit  finite 
difference  method.  For  the  simplification  of  the  calculation, 
the  diffusion  was  assumed  to  be  zero  at  the  inlet  and  the  out¬ 
let  of  the  channels  as  boundary  conditions.  This  condition 
could  be  experimentally  justified.  If  the  port  for  supplying 
the  gas  at  the  inlet  is  small  enough  to  have  high  speed,  and 
the  channels  are  narrow  enough  and  long  enough  without 
electrodes  at  the  downstream  of  the  cell  reaction  regions, 
the  resulting  convection  molar  flux  is  much  greater  than  the 
diffusion  molar  flux.  The  stationary  solution  for  the  direct 
current  case  was  calculated  by  assuming  the  open-circuit  as 
an  initial  condition.  The  result  was  then  used  as  an  initial 
condition  for  the  case  where  the  ac  current  superimposed 
on  the  dc  current.  The  load  current  /l  was  calculated  by 
integrating  the  output  current  density  Jo  over  the  entire 
cell  area  as  a  function  of  time.  The  Fourier  transforms  of 
the  load  current  and  the  applied  ac  voltage  were  used  to 
calculate  the  total  cell  impedance. 

The  calculation  conditions  of  the  reference  case  and  the 
parametric  study  are  shown  in  Table  1.  These  conditions 
were  based  on  the  experimental  setup  for  the  impedance 
measurements  at  Mitsubishi  Materials  Corporation  [6].  Ad¬ 
ditional  parameters  are  dc  voltage,  load  current  and  fuel  uti¬ 
lization.  Because  these  parameters  are  not  independent,  the 
fuel  utilization  factor  Uf  defined  below  will  be  mainly  used 
as  follows: 


(10) 


where  /l  the  load  current,  mu2  molar  flow  rate  of  the  sup¬ 
plied  hydrogen.  The  calculation  was  carried  out  with  an  ap¬ 
plied  ac  voltage  of  1  mV  amplitude  for  a  range  of  frequen¬ 
cies  between  100  and  0.01  Hz.  The  cell  was  meshed  with 
100  elements,  and  time  step  At  was  chosen  between  0.1  and 
5  ms  depending  on  the  frequency.  The  diffusion  coefficient 
was  evaluated  by  the  Chapman  and  Enskog  equation  as  a 
function  of  temperature  and  pressure  [13]. 


Table  1 

Simulation  conditions 


Parameter 

Reference 

case 

Parametric  study 

Cell  diameter  Dceu  (cm) 

8.5 

17 

Channel  effective  height  (cm) 

hf 

0.1 

0.05 

K 

0.1 

0.05 

Fuel  flow  rate  Gf  (Ncc/min/cm2) 

3 

6 

Air  flow  rate  Ga  (Ncc/min/cm2) 

15 

30 

Temperature  T  (K) 

1023.15 

973.15 

Total  pressure  p  (MPa) 

0.1013 

0.2026 

Fuel  utilization  Uf  (%) 

80 

10-90 

Local  cell  impedance  Rce n  (Q  cm2) 

0.5 

Table  2 

Diffusion  coefficient  (cm2/s) 

Df 

7.06 

0,  6.48a,  3.53b 

Da 

1.65 

0,  1.52a,  0.825b 

Inlet  H2O  mole  fraction  in  fuel  (%) 

1 

1 

Inlet  H2  mole  fraction  in  fuel  (%) 

99 

99 

Inlet  O2  mole  fraction  in  air  (%) 

20.95 

20.95 

Inlet  N2  mole  fraction  in  air  (%) 

79.05 

79.05 

a  For  973.15  K  at  0.1013  MPa. 
b  For  1023.15  K  at  0.2026  MPa. 


Here,  the  Peclet  number,  Pe  =  vL/D  ( v :  velocity,  L  the 
characteristic  length;  D  the  diffusion  coefficient),  which  is 
defined  as  the  ratio  of  the  convection  flux  and  the  diffusion 
flux,  is  estimated  to  be  about  2.4  for  the  radial  fuel  flow  and 
about  50  for  the  radial  air  flow,  by  taking  the  cell  outer  radius 
r0  (r0  =  4.25  cm)  as  the  characteristic  length  and  the  outlet 
velocity  vY0  (vV0  ~  4  cm/s  for  fuel  flow  and  vYO  ~  20  cm/s 
for  air  flow  at  the  reference  case)  as  the  velocity.  This  means 
that  the  convections  are  dominant  in  the  r-direction,  and  the 
r-direction  diffusion  in  the  fuel  channel  is  not  negligible. 
For  the  Peclet  number  in  the  z-direction,  Pe  =  0  for  the 
fuel  flow  channel  because  there  is  no  total  molar  flow  in 
the  z-direction  in  the  fuel  channel,  and  Pe  <  0.013  for  the 
air  flow  channel  with  the  cell  reaction  of  1  A/cm2  ( vz  = 
0.22  cm/s,  which  is  caused  by  the  consumption  of  O2  at 
the  electrode  surface).  That  is,  the  diffusion  is  much  larger 
than  the  convection  in  the  z-direction.  Further  more,  the 
parameter  93  =  vmh2 /r0D  is  estimated  to  be  93  ~  1.3  x 
10~3  for  the  fuel  flow  channel  and  93  ~  2.9  x  10-2  for 
the  air  flow  channel.  This  means  that  the  lateral  change 
in  the  mole  fraction  is  much  smaller  than  the  longitudinal 
change  at  steady  state.  The  parameter  Td  for  the  diffusion 
characteristic  time  in  the  z  direction  is  estimated  to  be  Td  = 
1.4  x  10-3  s  for  the  fuel  flow  channel  and  Td  =  6.1  x  10-3  s 
for  the  air  flow  channel.  Thus,  these  estimations  support  the 
assumption  that  the  perfect  mixing  in  the  z-direction  is  valid 
for  the  impedance  simulation  at  the  frequency  below  26  Hz 
(l/(27ZTd)).  If  the  impedance  due  to  gas  diffusion  in  the 
z-direction  at  frequencies  higher  than  1  /(27TTd)  is  important, 
such  impedance  can  be  considered  as  a  component  of  the 
local  impedance  as  well  as  the  impedances  for  the  finite  rate 
reaction,  etc. 
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Fig.  3.  Radial  distributions  of  electromotive  force,  current  density,  and 
partial  pressures  of  H2,  H2O  and  O2  under  the  reference  conditions.  Solid 
line:  with  diffusion,  broken  line:  without  diffusion. 


3.  Results  and  discussion 

3.1.  Voltage-current  characteristics  and  differential 
resistance 

The  steady-state  distributions  of  the  electromotive  force, 
current  density  and  partial  pressures  of  chemical  species  for 
the  fuel  utilization  of  80%  under  the  reference  conditions  are 
shown  in  Fig.  3.  Due  to  the  current  flow  dictated  by  the  cell 
reaction,  partial  pressure  values  of  the  hydrogen  and  oxygen 
decrease,  while  the  steam  partial  pressure  increases  with  the 
gas  flow.  As  a  result,  the  electromotive  force  and  the  current 
density  decrease  in  the  radial  direction.  The  results  obtained 
by  neglecting  the  diffusion  under  the  reference  conditions 
are  also  shown  in  Fig.  3.  By  comparing  the  two  cases,  it  is 
noticeable  that  the  increase  in  water  vapor  along  the  flow 
direction  is  promoted  by  the  diffusion  making  the  water  par¬ 
tial  pressure  higher  and  the  hydrogen  partial  pressure  lower 
in  the  upper-stream  region.  Thus,  the  drop  in  the  electromo¬ 
tive  force  becomes  larger  along  the  flow  direction,  and  the 
cell  voltage  decreases  inevitably  in  order  to  output  the  de¬ 
manded  value  of  the  total  current  (fuel  utilization).  In  other 
words,  the  cell  performance  gets  lowered  by  the  existence 
of  diffusion. 

Under  the  reference  conditions,  the  relationships  between 
the  cell  voltage  and  the  fuel  utilization  (or  the  load  cur¬ 
rent  density)  and  the  derivative  of  the  cell  voltage  with  re¬ 
spect  to  the  current  density  are  shown  in  Fig.  4.  The  cell 
voltage  decreases  with  current  density,  and  the  change  in 
the  cell  voltage  becomes  greater  as  the  fuel  utilization  ap¬ 
proaches  zero  or  100%.  This  occurs  because  the  depen¬ 
dence  of  the  electromotive  force  to  partial  pressures  becomes 
greater  as  the  fuel  utilization  approaches  zero  or  100%. 
The  derivative  of  the  cell  voltage  with  respect  to  the  cur¬ 
rent  density  has  the  dimension  of  resistance  per  unit  area, 
and  it  will  be  henceforth  referred  to  as  the  differential  cell 
resistance. 
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Fig.  4.  Variations  of  cell  voltage  Vceii  differential  cell  resistance 
dVceii/d/ave  and  real  parts  of  calculated  total  cell  impedance  at  100  and 
0.01  Hz  with  fuel  utilization  for  the  reference  conditions. 

3.2.  impedance 

The  variations  of  the  calculated  whole-cell  impedance 
versus  the  fuel  utilization  for  frequencies  between  100 
and  0.01  Hz  under  the  reference  conditions  are  shown  as 
Cole-Cole  plots  in  Fig.  5.  It  is  very  interesting  to  note  that 
the  Cole-Cole  plots  depict  perfect  semicircles  despite  the 
assumption  of  a  pure  resistance  as  local  impedance  of  the 
cell.  This  describes  the  change  in  electromotive  force  due  to 
the  cell  current  through  the  change  in  the  reactant  concen¬ 
tration  by  the  cell  reaction.  The  impedance  with  the  same 
mechanism  in  an  ideal  cell  with  a  small  volume  has  already 
been  reported  as  gas  conversion  impedance  by  Primdahl 
and  Mogensen  [14].  By  assuming  a  perfect  mixing  in  the 
gas  related  to  the  anode  reaction,  they  analyzed  the  change 
in  the  electromotive  force  caused  by  the  gas  conversion  due 
to  the  cell  reaction  in  the  anode  and  analytically  derived 
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Fig.  5.  Impedance  Cole-Cole  plots  between  0.01  and  100  Hz  with  param¬ 
eter  of  fuel  utilization  under  the  reference  conditions.  Open  and  closed 
symbols:  simulation  results;  Solid  line:  fitting  by  an  equivalent  circuit 
with  an  RC  parallel  circuit  as  shown  in  Fig.  7a.  Closed  symbol  shows 
the  point  of  the  frequency  value  noted  in  the  figure. 


K.  Takano  et  al.  /  Journal  of  Power  Sources  132  (2004)  42-51 


47 


an  expression  for  the  gas  conversion  impedance  as  an  RC 
parallel  circuit.  The  ideal  cell  with  perfect  mixing  should 
be  sufficiently  small  with  enough  diffusion,  and  it  is  just  an 
entrance  element  of  the  segmented  cell  in  our  simulation. 
The  impedance  of  the  whole  cell  cannot  be  obtained  by  a 
simple  integration  over  the  grid.  Therefore,  the  simulation 
is  an  effective  method  to  obtain  the  total  cell  impedance. 

For  the  simulation  results  shown  in  Fig.  5,  the  high- 
frequency  edge  of  the  semicircle  does  not  change  even  by 
changing  the  fuel  utilization.  However,  the  diameter  of  the 
semicircle  changes  with  the  fuel  utilization,  and  it  becomes 
bigger  as  the  fuel  utilization  approaches  either  zero  or  100%, 
having  a  minimum  value  at  the  fuel  utilization  of  about 
50-60%. 

Such  impedance  characteristics  can  be  explained  as  fol¬ 
lows:  When  the  applied  alternating  current  has  a  frequency 
enough  higher  than  the  transit  time  of  the  gas  through  the 
cell,  the  concentration-change  of  the  chemical  species  by  the 
applied  alternating  current  becomes  nearly  zero,  because  it 
is  proportional  to  the  integrated  value  of  the  reaction  quan¬ 
tity  in  a  half  cycle.  Any  change  in  the  electromotive  force 
does  not  occur  if  the  partial  pressures  of  the  chemical  species 
remain  constant.  Thus,  the  total  cell  impedance  agrees  with 
the  integrated  local  impedance.  On  the  other  hand,  for  the 
applied  alternating  current  with  low  frequency  which  can 
be  considered  as  a  small  change  in  the  direct  current,  the 
impedance  agrees  with  the  previously  mentioned  differen¬ 
tial  resistance  of  dc,  V-I  characteristics.  This  fact  is  indi¬ 
cated  as  the  simulation  result  in  Fig.  4,  in  which  the  real 
parts  of  the  impedances  at  the  high-frequency  edge  and  at 
the  low-frequency  edge  are  shown,  as  well  as  the  local  cell 
resistance  and  the  differential  resistance.  Therefore,  the  size 
of  the  semicircle  changes  with  the  fuel  utilization  according 
to  the  change  in  the  differential  cell  resistance. 

Fig.  5  also  shows  the  fitted  curves  superimposed  onto 
the  simulation  results  using  the  RC  parallel  circuit.  It  is 
very  interesting  to  notice  that  the  whole-cell  gas  conversion 
impedance  of  a  large-size  cell  can  also  be  successfully  ex¬ 
pressed  by  an  RC  parallel  circuit  used  for  a  small- size  ideal 
cell. 

According  to  Primdahl  and  Mogensen  [14],  the  frequency 
(that  will  be  henceforth  referred  to  as  the  characteristic  fre¬ 
quency)  for  which  the  imaginary  part  of  impedance  has  a 
maximum  value,  /g  is  given  by 


N 

2rtV 


(11) 


where  N  the  fuel  volume  flow  rate,  V  the  separate  volume 
related  with  the  anode  reaction.  The  equation  dictates  that 
the  characteristic  frequency  is  determined  from  the  period 
in  which  the  volume  V  of  the  gas  is  replaced  by  the  flow, 
that  is,  the  transit  time  of  gas  through  the  cell.  Regarding 
the  reference  conditions,  the  period  calculated  by  using  the 
volume  of  the  fuel  channel  as  V  is  0.53  s,  and  the  character¬ 
istic  frequency  is  0.3  Hz.  This  value  is  about  the  same  order 
of  magnitude  but  lower  than  the  characteristic  frequency  in 


the  simulation  results  depicted  in  Fig.  5.  The  reason  for  this 
seems  to  be  as  follows:  There  is  a  similar  phenomenon  in  the 
cathode  reaction  as  the  gas  conversion  in  the  anode  reaction, 
and  the  characteristic  frequency  for  the  cathode  is  estimated 
to  be  1.5  Hz  for  the  air  flow  channel.  The  simulation  results 
in  Fig.  5  are  convolutions  of  two  phenomena  for  the  anode 
and  cathode.  Here  it  should  be  remarked  that  the  influence 
on  the  impedance  of  the  cathode  phenomenon  is  relatively 
small  because  the  change  in  the  oxygen  partial  pressure  in 
the  air  flow  channel  due  to  the  cell  current  is  small,  and  the 
electromotive  force  change  is  small.  Fig.  6  shows  the  sim¬ 
ulation  result  for  the  case  without  diffusion  under  the  refer¬ 
ence  conditions.  This  indicates  that  the  diffusion  in  the  radial 
direction  causes  the  characteristic  frequency  to  be  lower. 

The  expressions  of  the  gas  conversion  impedance  given 
by  Primdahl  and  Mogensen  [14]  are 


RT\2  /  1  1  \  1 

2 F)  Vph2o  Ph2)  N 
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where,  Rg  and  Cg  are  the  resistance  and  the  capacitance  of 
a  parallel  circuit  for  the  gas  conversion  impedance,  respec¬ 
tively.  These  equations  are  valid  for  the  qualitative  discus¬ 
sions  of  the  effect  of  various  parameters  on  the  impedance 
behaviors  of  the  cells  in  large  scale  as  will  be  mentioned 
in  the  following  section.  However,  in  the  practical  loaded 
conditions,  they  cannot  be  applied  for  the  evaluations  of 
the  impedance  of  the  large  cells  where  the  concentrations 
of  gaseous  species  are  distributed  and  have  no  specific 
values. 
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Fig.  6.  Impedance  Cole-Cole  plots  between  0.01  and  100  Hz  with  pa¬ 
rameter  of  fuel  utilization  for  the  reference  conditions  without  diffusion. 
Open  and  closed  symbols:  simulation  results;  Solid  line:  fitting  by  an 
equivalent  circuit  with  an  RC  parallel  circuit  as  shown  in  Fig.  7a.  Closed 
symbol  shows  the  point  of  the  frequency  value  noted  in  the  figure. 
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Fig.  7.  Equivalent  circuit  for  the  whole-cell  impedance  considering  the 
gas-conversion  impedance  by  an  RC  parallel  circuit  with  Rwg  and  Cwg. 

3.3.  Influence  of  design  parameters  on  total  cell 
impedance 

The  influence  of  various  design  parameters,  such  as 
cell  dimensions  and  operating  conditions,  on  the  total  cell 
impedance  was  parametrically  examined.  Parameter  fitting 
was  performed  for  the  simulation  results  by  the  impedance 
equivalent  circuit  depicted  in  Fig.  7.  This  equivalent  cir¬ 
cuit  for  the  alternating  current  was  derived  from  the  SOFC 
equivalent  circuit  of  Fig.  2  by  substituting  an  RC  parallel 
circuit  with  Rwg  and  Cwg  for  the  electromotive  force.  In 
addition,  the  total  cell  impedance  was  normalized  by  the 
cell  area.  The  error  in  each  fitting  was  sufficiently  small 
as  well  as  in  Figs.  5  and  6,  and  it  was  confirmed  that  the 
impedance  on  the  high-frequency  edge  agreed  with  the 
local  cell  resistance  in  each  case. 

Figs.  8  and  9  show  Rwg  and  the  time  constant  of  rwg  = 
RWgCWg  for  the  whole-cell  gas-conversion  impedance  for 
the  reference  case  and  for  the  cases,  such  as  neglected  diffu¬ 
sion,  double  gas  flow  rate,  half  channel  height,  double  cell 
diameter,  temperatures  lowered  to  973  K,  or  pressure  dou¬ 
bled  under  the  reference  conditions. 

The  influence  of  most  of  the  parameters  on  the  whole-cell 
gas-conversion  impedance  can  be  explained  in  principle  via 
Eqs.  (11)  and  (12). 
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Fig.  8.  Rw g  vs.  fuel  utilization  with  changing  parameter:  (a)  reference  case; 
(b)  neglected  diffusion;  (c)  doubled  gas  flow  rates;  (d)  halved  channel 
heights;  (e)  doubled  cell  diameter;  (f)  973  K  temperature  and  (g)  doubled 
total  pressure. 


Fig.  9.  rwg  vs.  fuel  utilization  with  changing  parameter:  (a)  reference  case; 
(b)  neglected  diffusion;  (c)  doubled  gas  flow  rates;  (d)  halved  channel 
heights;  (e)  doubled  cell  diameter;  (f)  973  K  temperature  and  (g)  doubled 
total  pressure. 

For  all  cases  in  Fig.  8,  normalized  resistance  Rwg  is  a 
function  of  the  fuel  utilization  and  becomes  bigger  as  the 
fuel  utilization  approaches  either  0  or  100%,  having  a  mini¬ 
mum  value  at  the  fuel  utilization  around  50-60%.  This  may 
be  guessed  from  Eq.  (12)  by  considering  the  dependence 
of  Rg  on  the  partial  pressures  of  hydrogen  and  water  vapor. 
As  shown  in  Fig.  9,  the  time  constant  rwg  scarcely  changes 
with  the  fuel  utilization  between  20  and  80%,  but  consid¬ 
erably  increases  as  the  fuel  utilization  approaches  zero  or 
100%,  although  the  transit  time  of  the  gas  does  not  change. 
Such  behavior  cannot  be  simply  explained,  but  it  qualita¬ 
tively  agrees  with  the  experimental  results  [6] .  It  seems  to  be 
related  to  the  distribution  of  the  chemical  species  in  the  cell. 

When  the  diffusion  is  neglected,  the  cell  performance  is 
improved,  resulting  in  a  reduced  Rwg  with  a  decrease  in  time 
constant  as  previously  mentioned.  When  the  gas  flow  rate 
is  doubled  or  when  the  channel  height  is  reduced  to  a  half 
of  the  current  value,  then  rwg  decreases  less  than  a  half, 
although  the  transit  time  of  the  gas  becomes  50%  shorter. 
This  seems  to  be  caused  by  a  relatively  weak  effect  of  dif¬ 
fusion,  as  the  mass  transfer  by  convection  doubles  due  to 
the  doubled  flow  velocity,  while  the  mass  transfer  by  dif¬ 
fusion  is  hardly  affected.  The  transit  time  does  not  change, 
even  if  the  cell  diameter  is  doubled  for  a  given  gas  flow 
rate  for  unit  cell  area.  However,  there  is  a  decrease  in  the 
time  constant,  which  can  also  be  explained  by  the  same  rea¬ 
son.  According  to  Eq.  (12),  the  gas  conversion  impedance 
is  inversely  proportional  to  the  gas  flow  rate;  however,  Rwg 
becomes  less  than  one-half  by  doubling  the  gas  flow  rate. 
This  over-decreasing  can  also  be  explained  by  the  relatively 
weak  effect  of  the  diffusion  by  doubling  the  gas  flow  veloc¬ 
ity.  For  the  same  reason,  Rwg  is  decreased  by  halving  the 
channel  height  or  by  doubling  the  cell  diameter. 

With  a  50  K  decrease  in  temperature  from  1023  K,  there 
appears  a  small  decrease  of  about  5-6%  in  Rwg  as  shown  by 
(f)  in  Fig.  8.  Most  part  of  the  decrease  can  be  explained  by  the 
fact  that  the  gas  conversion  impedance  is  proportional  to  the 
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Table  2 


Equivalent  circuit  parameters  for  the  parametric  study  of  Figs.  10-12 


Case 

a 

b 

c 

d 

e 

Equivalent  circuit  in  Fig.  2 

(a) 

(a) 

(b) 

(b) 

(b) 

Rceii  (£2  cm2) 

0.5 

1.0 

— 

— 

— 

Ri  cm2) 

— 

— 

0.5 

0.25 

0.25 

Rr  (12  cm2) 

— 

— 

0 

0 

0.1 

Cd  (F/cm2) 

— 

— 

0 

0 

0.16 

Re  (ficm2) 

— 

— 

12.5 

10.0 

10.0 

Rs  (12  cm2) 

— 

— 

0 

0.25 

0.15 

temperature,  as  dictated  by  Eq.  (12),  though  Rwg  is  affected 
by  a  decrease  in  the  diffusion  coefficients  with  decreased 
temperature.  While  the  time  constant  increases  about  4%  by 
the  decreased  temperature  as  shown  by  (f)  in  Fig.  9,  most 
of  this  increase  can  be  also  explained  by  Eq.  (11).  As  the 
total  pressure  increased  twice,  Rwg  does  not  change,  and  rwg 
is  doubled  as  shown  in  Figs.  8  and  9.  This  can  be  verified 
via  Eqs.  (11)  and  (12),  since  the  partial  pressures  become 
twice  by  doubling  the  pressure,  simultaneously  resulting  in 
the  halved  volume  flow  rate. 

3.4.  Influence  of  local  cell  equivalent  circuit  parameters 
on  total  cell  impedance 

To  examine  the  correlation  between  the  local  impedance 
and  the  total  cell  impedance,  a  simulation  has  been  car¬ 
ried  out  for  a  variety  of  equivalent  circuit  parameters  as 
shown  in  Table  2.  The  results  for  fuel  utilization  of  70  and 
80%  under  the  reference  conditions  are  shown  in  Figs.  10 
and  11. 

In  every  case  except  case  e,  the  normalized  impedance 
on  the  high-frequency  edge  agrees  with  the  local  cell 
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Fig.  10.  Cole-Cole  plots  at  70%  fuel  utilization  for  various  cases  (a,  b, 
c,  d,  e)  of  local-impedance  parameters  as  shown  in  Table  2.  Open  and 
closed  symbols:  simulation  results;  Solid  line:  fitting  by  an  equivalent 
circuit  with  an  RC  parallel  circuit  as  shown  in  Fig.  7b.  Closed  symbol 
shows  the  point  of  the  frequency  value  noted  in  the  figure. 
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Fig.  11.  Cole-Cole  plots  at  80%  fuel  utilization  for  various  cases  (a,  b, 
c,  d,  e)  of  local-impedance  parameters  as  shown  in  Table  2.  Open  and 
closed  symbols:  simulation  results;  Solid  line:  fitting  by  an  equivalent 
circuit  with  an  RC  parallel  circuit  as  shown  in  Fig.  7b.  Closed  symbol 
shows  the  point  of  the  frequency  value  noted  in  the  figure. 

impedance.  In  case  e,  considering  the  complex  impedance 
of  an  RC  parallel  circuit  with  the  characteristic  frequency 
of  10  Hz,  the  Cole-Cole  plot  looks  as  if  a  summation  of  the 
two  semicircle  impedances  with  characteristic  frequencies 
are  10  and  0.4  Hz.  In  this  case,  the  normalized  whole-cell 
impedance  on  the  high-frequency  edge  agrees  with  the  lo¬ 
cal  impedance  on  the  high-frequency  edge.  It  is  indicated 
that  the  local  impedance  directly  appears  as  the  total  cell 
impedance,  when  it  is  a  frequency-depending  complex 
impedance  at  much  higher  frequency  than  the  characteris¬ 
tic  frequency  of  the  gas-conversion  impedance.  It  was  also 
suggested  that  the  total  impedance  of  the  cell  appears  as  the 
convolution  of  the  local  impedance  and  the  gas  conversion 
impedance. 

When  the  local  cell  resistance  Rcen  increases  from 
0.5  ^  cm2  of  the  reference  value  to  1.0  £2  cm2,  the  semicir¬ 
cle  of  the  Cole-Cole  plot  only  parallel  moves  by  0.5  £2  cm2. 
However,  the  size  of  the  impedance  semicircle  increases 
in  all  cases,  when  the  electronic  conductance  is  intro¬ 
duced  using  the  equivalent  circuit  of  Fig.  2b  for  the  local 
impedance.  The  increase  in  the  impedance  is  greater  for 
80%  fuel  utilization  than  for  70%  utilization.  Such  increase 
in  the  impedance  may  be  related  to  the  increase  in  the  ef¬ 
fective  fuel  utilization  induced  by  the  current  leak  due  to 
the  electronic  conductance.  Therefore,  the  simulation  for 
the  effective  fuel  utilization  t/eff ,  which  is  defined  using 
the  total  reaction  current  instead  of  the  total  load  current  in 
Eq.  (10),  was  carried  out,  and  the  parameter  fitting  to  the 
simulation  results  using  the  equivalent  circuit  in  Fig.  7b  was 
performed.  As  a  result,  Rwg  and  rwg  are  shown  against  the 
effective  fuel  utilization  t/eff  in  Figs.  12  and  13.  Most  of  the 
plots  in  both  figures  are  overlapping  with  each  other,  which 
indicate  that  the  whole-cell  gas-conversion  impedance  (Rwg 
and  Twg)  for  a  fixed  effective  fuel  utilization  scarcely  de¬ 
pends  on  the  local  impedance. 
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Fig.  12.  Rw g  vs.  effective  fuel  utilization  for  the  various  local-impedance 
parameters  as  shown  in  Table  2  as  cases  a,  b,  c,  d  and  e. 


Fig.  13.  rwg  vs.  effective  fuel  utilization  for  the  various  local-impedance 
parameters  as  shown  in  Table  2  as  cases  a,  b,  c,  d  and  e. 


for  the  intermediate  values  of  the  fuel  utilization.  The 
whole-cell  gas-conversion  impedance  is  roughly  in¬ 
versely  proportional  to  the  supplied  gas  flow  rate,  being 
affected  by  the  diffusion  in  the  flow  direction. 

(3)  The  whole-cell  gas-conversion  impedance  can  be  well 
expressed  by  an  RC  parallel  circuit.  Its  time  constant  is 
roughly  equal  to  the  transit  time  of  the  fuel  gas  through 
the  cell,  being  affected  by  the  diffusion  in  the  flow  di¬ 
rection.  It  weakly  depends  on  the  fuel  utilization,  but 
considerably  increases  as  the  fuel  utilization  approaches 
either  zero  or  100% 

(4)  In  an  SOFC  cell  using  electrolyte  with  electronic  con¬ 
ductance,  it  is  necessary  to  evaluate  the  whole-cell 
gas-conversion  impedance  using  the  effective  fuel 
utilization  considering  the  current  leak  due  to  the 
electronic  conduction.  The  whole-cell  gas-conversion 
impedance  for  a  given  effective  fuel  utilization  weakly 
depends  on  the  local  cell  impedance. 

From  the  simulation  results  it  is  suggested  that  the  total 
impedance  of  the  cell  appears  as  the  convolution  of  the  local 
impedance  and  the  gas  conversion  impedance. 
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Appendix  A.  Discussion  on  the  parameter  *R 


4.  Conclusions 

A  numerical  simulation  method  which  links  the  local 
impedance  with  the  total  cell  impedance  at  power  genera¬ 
tion  has  been  proposed.  The  simulation  of  the  simple  cases 
under  experimental  conditions  is  carried  out  to  obtain  the 
total  cell  impedance  by  changing  parameters,  such  as  ac  fre¬ 
quency,  fuel  utilization,  gas  flow  rate,  etc.  The  following 
results  have  been  found: 

(1)  The  Cole-Cole  plot  for  the  total  cell  impedance  shows 
a  capacitive  semicircle,  even  if  a  pure  resistance  is  as¬ 
sumed  as  the  local  impedance,  due  to  the  gas  conver¬ 
sion  impedance  which  describes  the  change  in  the  elec¬ 
tromotive  force  by  the  kinetic  current  through  the  gas 
conversion  due  to  the  reaction.  The  high-frequency  edge 
in  the  Cole-Cole  plot  agrees  with  the  assumed  local 
impedance. 

(2)  The  whole-cell  gas-conversion  impedance  is  depen¬ 

dent  on  the  fuel  utilization,  supplied  gas  flow  rate,  cell 
dimension,  etc.  It  becomes  larger  as  the  fuel  utiliza¬ 
tion  approaches  zero  or  100%  and  has  minimum  value 


The  three-dimensional  molar-transport  equation  for  a 
species  M  when  no  spatial  species  production  is  given  as 
follows  [15]: 

~  +  V(cm?)  =  V(cDmVtm).  (A.l) 

where  cm  the  molar  density  (mol/cm3),  c  the  total  molar 
density  (mol/cm3),  v  the  molar  average  flow  velocity  vector 
(cm/s),  Dm  the  diffusion  coefficient  (cm1 2/s),  xm  the  mole 
fraction. 

When  axial  symmetry  and  constant  molar  density  are  as¬ 
sumed,  Eq.  (A.l)  can  be  described  as  follows: 


3vm 

~dT 


l  3 

+  -  —  (rvrxM)  + 

r  or 


dvz* M 
3  z 


1  3 

r  dr 


(A. 2) 


To  compare  the  distribution  of  the  mole  fraction  between 
that  in  the  z  direction  and  that  in  the  r  direction,  we  neglect 
the  r-direction  diffusion  term  and  the  z-direction  convection 
in  the  above  equation  and  assume  a  steady  state  for  simplic- 
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ity.  Using  continuity  equation,  the  previous  equation  can  be 
normalized  as  follows: 

3xm  d^X  M 

m  O=o 

(A3) 

where 

r 

$  =  — 
r0 

(A.4) 

Z 

? =h 

(A. 5) 

<n  = 

t o  Dm 

(A.  6) 

where  r0  is  the  outlet  radius  of  cell;  vro  the  radial  molar 
average  flow  velocity  at  cell  outlet;  h  the  height  of  flow 
channel. 

By  integrating  the  Eq.  (A. 3)  for  the  cell  volume,  the  fol¬ 
lowing  equation  is  obtained: 

/  (*M,£=1  —  *M,£=0) 

JO 

=  f  \  (xMf=\  —  YM,£=o)d£ df  (A.7) 

Jo  Jo 

where  is  the  mean  value  of  93. 

The  Eq.  (A.7)  indicates  that  the  difference  of  the  mole 
fraction  between  f  =  0  (electrode  surface)  and  f  =  1  (sepa¬ 
rator  surface)  is  93  times  the  difference  between  $  =  0  (inlet) 
and  $  =  1  (outlet)  as  the  average  values.  That  is,  the  lateral 
change  of  the  mole  fraction  in  the  channel  is  evaluated  to 
be  roughly  93  times  the  longitudinal  change  at  steady  state. 
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